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Abstract 

The Monte Carlo method in the canonical ensemble is used to investigate magnetization reversal 
in amorphous transition metal - rare earth multilayers. Our study is based on a model containing 
diluted clusters which exhibit an effective uniaxial anisotropy in competition with random mag- 
netic anisotropy in the matrix. We simulate hysteresis loops for an abrupt profile and a diffuse 
one obtained from atom probe tomography analyses. Our results evidence that the atom probe 
tomography profile favors perpendicular magnetic anisotropy in agreement with magnetic measure- 
ments. Moreover, the hysteresis loops calculated at several temperatures qualitatively agree with 
the experimental ones. 
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I. INTRODUCTION 



In the last decades, compositionally modulated films have been extensively investigated 
because of their interesting magnetic properties and potential applications in magneto- 
optical recording. In particular, transition metal (TM) / rare earth (RE) amorphous multi- 
layers, such as Fe/Dy multilayers, have been studied because of to their particular magnetic 
properties [3, S, li, 5, 6, 3,3,1^, 10] . In some conditions, they may exhibit a large magnetic 
anisotropy perpendicular to the layer, a high Curie temperature and a high coercive field for 
example. Nevertheless, the origin of the perpendicular magnetic anisotropy (PMA) is not 
yet clearly understood and can be explained by different models. These models are based 
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either on an anisotropic distribution of TM-RE pairs along the perpendicular direction 
dipolar interactions jlj], local structural anisotropy 3, [ijl or single ion anisotropy [7]. 
Recently, the results of Tamion et al. on (Fe 3nm/Dy 2nm) multilayers 0, [lo| have allowed 
to correlate the uniaxial anisotropy constant with the elaboration temperature by means 
of atom probe tomography, SQUID magnetometry and polarized neutron reflectivity mea- 
surements. It has been shown that the PMA is maximum for an elaboration temperature 
Ts = 570K, that is for diffuse interfaces. Consequently, PMA is rather due to the formation 
of an amorphous Fe-Dy alloy at the interfaces than abrupt interfaces. 

In the present study, we investigate the influence of the concentration profile along the 
multilayer on the PMA of amorphous Fe/Dy multilayers by means of Monte Carlo simu- 
lations of hysteresis loops. In the framework of a local structural anisotropy model [l^] 
we study the influence of the single-ion anisotropy constant on the RE atoms. Our aim is 
to reproduce qualitatively the experimental hysteresis loops. Thus, we are able to propose 
magnetization reversal mechanisms in relation with the local magnetic anisotropy and the 
concentration profile. Our numerical results also provide complementary information such 
as magnetization profile along the multilayers for different values of the applied field. The 
model and simulation technique are described in Sec. [Ill In the following sections, we report 
and discuss numerical results on the influence of the concentration profile at low temperature 
for different types of magnetic anisotropy: in the case of uniaxial anisotropy on all Dy sites 
in Sec. IIH] for the cluster anisotropy model (described below) in Sec. [TV] and in the case of 
the cluster anisotropy model combined with random anisotropy in Sec. |VJ More precisely, 



Sec. IIHI is devoted to a simple model for which theoretical results are available which has 
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allowed us to validate the simulation technique. The investigation, in Sec. HVl of the cluster 
anisotropy model is performed to understand the differences with the uniaxial anisotropy 
model. The results of this section are very useful to explain the numerical data obtained in 
Sec. El In Sec. IVI[ we report the temperature effect on hysteresis loops and compare our 
results to experimental ones. Finally a conclusion is given in Sec. IVIII 



II. MODEL AND MONTE CARLO SIMULATION 



A. Description of the model 

Our model consists in a face centered cubic (FCC) multilayer system made up of Fe and 
Dy atoms. We choose this closed-packed structure (each atom has 12 nearest neighbors) 
because its density is very close to those of amorphous structures 15, Q- To take into 
account the different interatomic distances observed in an amorphous material, we consider 
distributed exchange interactions since they strongly depend on the interatomic distances 

Q. 

It has been shown previously that atomic diffusion has a majour influence on the macro- 
scopic magnetic anisotropy [9]. Consequently, two different concentration profiles are used 
to explain these features. The first one, called abrupt in the following, corresponds to a 
multilayer made up of pure Fe and Dy layers with abrupt interfaces (Fig. [TJ(a)). The second 
profile has been directly obtained from atom probe tomography analysis of (Fe 3nm/Dy 
2nm) multilayers elaborated at 570K it is called APT (atom probe tomography) profile 
in the following. This profile is composed of a Fe-rich region (FegoDyio) and a region in 
which the concentration varies; we would like to note that this profile does not display any 
pure region (Fig. (U(b)). As it has been previously mentioned, this concentration profile 
leads experimentally to the maximum of PMA. In our model, there are 36 atomic planes 
corresponding to a double (Fe 3nm/Dy 2nm) bilayer. For simplicity, we do not take into 
account the difference in atomic radius of Fe and Dy atoms. We apply periodic boundary 
conditions in the film plane and free boundary conditions in the ^-direction. In sections |V] 
and IVIl in order to compare with experimental results, we calculate the physical quantities 
over the central bilayer only, to get rid off free surface effects. 

Each node of the FCC lattice is occupied by a classical Heisenberg spin which is a 3D 
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FIG. 1: Concentration profile of each specie along the multilayer growth direction for an abrupt 
(a) and a APT (b) concentration profile. 



vector spin. This classical model is the most suitable for investigating magnetic configu- 
rations in systems with competitive interactions (exchange, anisotropy and Zeeman). We 
consider the following Hamiltonian: 



E = ~ £ ^(Si.S^ + K-B^nii, (1) 

<i,j> i 

where Sj is a classical Heisenberg spin and nij is the magnetic moment at site i. The first term 
corresponds to exchange energy. Here, we consider nearest-neighbor (NN) interactions only, 
so Jij are the NN exchange interactions. The second term, E & , is the magnetic anisotropy 
energy. We consider only single-ion anisotropy on the Dy atoms, the magnetic anisotropy 
of the Fe atoms is neglected. The third term describes the Zeeman energy where B is the 

n 

applied magnetic field. For the Fe spins, we take the values proposed by Heiman et al. [16| 
which depend on the local concentration under the form: 

S Fc (X Fc ) = 1.1- 1.125(1 - X Fc ) (X Fc > 0.4). (2) 

Since the Dy spin is not sensitive to the local environment, we take Sby = 2.5. The Lande 
factors are g Fe = 2 and go y = 4/3. The atomic moments for Fe and Dy are related to the 
spins according to: 



m Fe — gFe/^BSW^Fe) 
9By 



9D> 



I^bSdy — 10/iB- 



(3) 
(4) 



The NN exchange interactions are also extracted from the work of Heiman et al. [if]] on 
Fe-Dy amorphous thin films. Since they have been obtained by mean-field calculations, we 
have to adjust them to get the pure amorphous Fe and Dy Curie temperatures (respectively 
270K (aQ and 110K [5]) by Monte Carlo simulations. 

The Fe-Fe and Fe-Dy exchange interactions linearly depend on the local concentration 
Xp e . The Dy-Dy exchange interaction, which is much smaller than the others, is independent 
of the concentration. These exchange interactions are given by: 

J Fe _ Fc (X Fe )/A;B = 77 + 449(1 - X Fc ) (in K), (5) 

J Fc _ Dy (X Fe )/A;B = 8 - 198(1 - X Fc ) (in K), (6) 

^Dy-Dy/^B = 6.5 (in K). (7) 

J Fe _D y is strongly negative; it is thus responsible for the ferrimagnetic order (with Fe mo- 
ments antiparallel to Dy moments) which has been experimentally observed by polarized 
neutron reflectivity measurements [10| . 

In the following, we define the ferrimagnetic magnetization per atom (Mf err i magne tic) by: 

-Mferrimagnetic = "T7" /",(! — -^ F e) m Dy ~ X Fe m Fe , (8) 
iV P 1=1 

where N p is the number of atomic planes and m Fe is the Fe atomic moment in the I th plane. 
For the abrupt and APT profile, Mf errimagnetic is equal to 1.87/xb and 1.21/i B respectively. 



B. Monte Carlo simulation 



In this work, we use Monte Carlo (MC) method in the canonical ensemble based on the 
Metropolis algorithm [19J. This famous method is known to ensure a good convergence 
towards the thermodynamic equilibrium but it may in some circumstances not describe 
properly the magnetization reversal. So we use here the algorithm proposed by Nowak et 
al. (2^]. In this algorithm, a site i is chosen randomly and the new orientation of the spin 
of this site within a cone with a given size is proposed. To achieve this procedure, a random 
vector u with uniform probability distribution within a sphere of radius R< 1 is added to 
the initial normalized spin Sj = Sj/ || Sj ||. The new orientation of the spin is then given by 
the unit vector: 
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The energy difference AE between the new (SJ and initial (Si) spin orientations is 
calculated from Eq.fJTJ (with S- =|| Sj || s- ). If AE < 0, the transition is accepted 
since it lowers the energy. Otherwise, the transition is accepted with a probability rate 
proportional to exp(— AE/k^T). One Monte Carlo step (MCS) consists in examining all 
spins of the system once. The variation of R allows varying the acceptance rate in order 
to optimize the efficiency of the algorithm. With this technique, reliable MC simulations 
with a reasonable number of MCS can be performed at low temperature thanks to large 
enough acceptance rates. Moreover, unlike the standard Metropolis algorithm, this method 
prevents non physical spin flips by tunneling accross the barrier. 

The initial magnetic configuration is completely disordered. Then, a strong external 
field +B max is applied along a given axis. The hysteresis loop is performed by varying the 
external field with a step AB = 0.1 5T. At each field, 40 000 MCS are discarded for local 
equilibration before averaging the magnetization components over the next 40 000 MCS. As 
we are interested in multilayers with configurational disorder, the final hysteresis loops are 
measured by averaging over several disorder configurations. 

In order to choose reasonable values of R, we perform numerical simulations on simple 

n 

systems which satisfy the Stoner-Wohlfarth (SW) model [21[. Then, we keep for R the values 
which allow to reproduce the theoretical hysteresis loops within a reasonable computational 
time (R ~ 0.1). 

III. UNIAXIAL ANISOTROPY 

In this section, we consider uniaxial anisotropy perpendicular to the film plane for all Dy 
atoms. The anisotropy energy is given by: 

E a = -D Dy Yl ( S - z ) 2 > ( 10 ) 

i£Dy sublattice 

where Dd y is the single- ion anisotropy constant. We investigate here the influence of the 
concentration profile and of the anisotropy constant on the magnetization reversal of an 
amorphous Fe/Dy multilayer. 
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FIG. 2: (a) Hysteresis loops in perpendicular and parallel orientations for an Fe/Dy multilayer with 
the abrupt profile and uniaxial magnetic anisotropy at T = IK for D-Qy/k-Q = IK. (b) Coercive 
field as a function of the Dy uniaxial anisotropy constant at T = IK. 

A. Abrupt profile 



The hysteresis loops of the multilayer with an abrupt profile at low temperature for an 
applied field in the film plane and perpendicular to it are shown in Fig. [21(a) . For the latter 
case, the hysteresis loop reveals two types of behaviors. For \B\ < Bq = 2,45 ± 0, 15T, the 
magnetic configuration is ferrimagnetic (Fe and Dy layers are antiparallel to each other). 
Beyond Bq, the absolute value of the magnetization increases with the external field. These 
configurations exhibit an interface domain wall (IDW) located in the Fe layer which has a 
free surface. This is due to the competition between the external field and the antiparallel 
Fe-Dy coupling. It has to be noted that an applied field of 6T is not strong enough to 
saturate the magnetization of the sample. For an in-plane field, the loop is closed and 
the absolute value of the magnetization increases monotonously with the field intensity as 
expected for a hysteresis loop measured in a hard orientation. 

The switching field predicted by the SW model, i.e. for uniform rotation process, is: 

? sw _ 2X Dy L) Dy S , D y 



M, 



(11) 



ferrimagnetic 



where Xd 7 is the Dy atomic fraction. In Fig. [2j(b), we have plotted the coercive field as 
a function of the Dy anisotropy constant for a perpendicular applied field. The simulated 
coercivity increases as expected with -Doy, but it is smaller than the SW values except 
for Dvy/k-Q < 0.25K. This deviation from the SW model confirms that the magnetization 
reversal is not uniform. The reason is that the antiferromagnetic Fe-Dy interactions favor 
the formation of the IDW during the reversal. Moreover, we note the existence of a spin-flop 
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FIG. 3: (a) Hysteresis loops in perpendicular and parallel orientations for an Fe/Dy multilayer 
with the APT profile and uniaxial magnetic anisotropy at T = IK for D^y/k-Q = IK. (b) Coercive 
field as a function of the Dy uniaxial anisotropy constant at T = IK. 

transition at 5sf ~ 2, IT for .Drjy/^B > IK in good agreement with previous studies on 
ferrimagnetic multilayers 22|, |23|]- The spin- flop transition corresponds to the formation of 
a twisted spin structure at the interface. 



B. APT profile 

In the case of the APT profile, the hysteresis loops at T = IK are shown in Fig. [31(a). 
Contrary to the previous case, we observe a square loop when the field is normal to the 
layers which characterizes a ferrimagnetic order for each value of the magnetic field. These 
results show that for the APT profile the system is equivalent to a homogeneous ferrimag- 
netic material (at least in the field range investigated). In Fig. [31(b), we have plotted the 
coercivity at T = IK as a function of the anisotropy constant -Do y - We obtain a reason- 
able agreement with the SW model up to D-£, y /k-B ~ IK, i.e. the magnetization reversal is 
roughly uniform. For D-^j/k-Q > IK, our results deviate significantly from the SW model. 
Indeed, the magnetization reversal starts in Fe-rich planes by breaking Fe-Dy antiparallel 
couplings which leads to a non-uniform rotation process. Then, we can conclude that the 
APT profile favors uniform rotation process in comparison with the abrupt profile. 
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FIG. 4: Thermal variation of the reduced magnetization (continuous line) and of the coercive field 
(dashed line) for D^ y /kB = IK ((a) abrupt profile; (b) APT profile). 

C. Influence of the temperature 

Finally, we have studied the temperature influence on the hysteresis loops for each con- 
centration profile. The thermal variation of the reduced magnetization and of the coercive 
field is shown in Fig. HI 

Our results evidence a faster decrease of the coercivity as the temperature increases in 
the case of the abrupt profile (Fig. Hj). The temperature dependence of the coercivity can be 
fitted to a square law contrary to the APT profile for which the decrease is roughly linear. 
This faster decrease is related to the decrease of the magnetization which is also faster than 



for the APT profile [24]. 



IV. CLUSTER ANISOTROPY 



Here, we consider the local anisotropy model. This model is based on the existence 
of small crystallized clusters at the scale of a few interatomic distances depending on the 
elaboration method. These clusters define on average an easy axis perpendicular to the film 
plane [3]. In our model, they consist of 13 atoms (one Fe central atom and its 12 nearest 
neighbors) (Fig. [5]). Among the 8 neighbors which are not in the xy plane of the central 
atom, between 2 and 4 Dy atoms are randomly distributed in order to obtain a cluster 
composition close to that of defined compounds. The easy axes (on Dy atoms belonging to 
the clusters) are along the Fe-Dy bonds; so the unit vectors z» of these 4 axes are (±^g,0,^=) 
or (0, ±-4=,4=). This leads on average to an easy axis perpendicular to the film plane. 
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FIG. 5: Schematic representation of the ordered clusters. Fe atoms are in blue, Dy atoms in green. 
The arrows show the local anisotropy axes on the Dy sites. 



In this case, the anisotropy energy is defined as: 



E, = -D Dy Yl ( S - z *) 2 - ( 12 ) 

i£Dy in clusters 

As the anisotropy coefficient D-q 7 in amorphous multilayers is not accurately determined in 
the literature, we consider it as a free parameter. In the same way, the cluster concentration 
(cciuster), which is defined as the number of atoms (Dy and Fe) included in the clusters 
divided by the total number of atoms, is a free parameter in the simulations since it cannot 
be evaluated experimentally. 

We investigate the influence of the Dy anisotropy constant and of the cluster concentration 
for the two concentration profiles. The cluster concentration varies from 5% to 20%, that 
is only 3.4% to 13.8% of the Dy atoms display single-ion anisotropy. Consequently, we 
take larger values of the anisotropy constant than in the previous section in order to obtain 
significant values of the coercive field. In the framework of the SW model, it can easily be 
shown that the coercive field is given by: 

^cluster j~\ c2 

B sw = *b, OpA» (i3) 

■^-'"ferrimagnetic 

where X^ ter is the atomic fraction of Dy atoms belonging to the clusters. 



A. Abrupt profile 

For the abrupt profile, the cluster concentration is fixed at 5%; it cannot go further 
because the clusters are located at the interfaces. Fig. [6j(a) shows the hysteresis loops 
for a magnetic field perpendicular to the film plane at T = IK and different values of 
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FIG. 6: (a) Hysteresis loops in the perpendicular orientation for an Fe/Dy multilayer with an 
abrupt profile at T = IK and different values of the Dy anisotropy constant ; (b) Simulated and 
SW coercive field variation as a function of the Dy anisotropy constant value. 

D-Dy/k-Q- We observe that the hysteresis loops are almost square indicating that the normal 
orientation is an easy axis. For all anisotropy constant values, the magnetic configuration is 
ferrimagnetic with the Dy moments being aligned with the external field. Unlike the case of 
uniaxial anisotropy examined in the previous section, we do not observe here IDW because 
the coercive field is always smaller than the spin-flop field (~ 2. IT). 

Fig. 0(b) shows the coercive field variation as a function of the Dy anisotropy constant 
value. Our results are in good agreement with the SW model up to D-£> y /k-Q = 30K which 
evidences a magnetization reversal by uniform rotation. Beyond this value, the reversal is 
no more purely uniform due to the clusters. 

B. APT profile 

In the case of the APT profile, the clusters are not only localized at the interfaces. The 
hysteresis loops are represented on Fig. [3(a) at T = IK for D^y/k-Q = 10K and several 
cluster concentration values. As for the abrupt profile, our results indicate that the z-axis 
is an easy orientation and no IDW is observed for the same reason. 

As can be seen in Fig. 0(b), for D^y/k-Q = 10K, the magnetization reversal process 
corresponds to an uniform rotation for all values of c c i uste r, whereas for D Dy //c B = 50K the 
magnetization reversal is not purely uniform anymore starting from c c i us t e r ~ 15%. 

The coercive field is proportional to the cluster concentration and to the Dy anisotropy 
constant as far as the product c c i uster .D Dy //cB remains small enough (< 5K). We would like 
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FIG. 7: (a) Hysteresis loops in the perpendicular orientation for an Fe/Dy multilayer with an 
APT profile at T = IK and for D^ y /k-Q = 10K and different cluster concentrations; (b) Simulated 
and SW coercive field variation as a function of the cluster concentration for D-£, y /k-Q = 10K and 
D Dy /k B = 50K. 

to note that, for a given anisotropy constant, the coercive field is slightly larger for an APT 
profile than for an abrupt profile (Tab. HI). 

V. CLUSTER AND RANDOM MAGNETIC ANISOTROPY 

In this section, we study a more realistic model of Fe/Dy multilayers to propose a qual- 
itative explanation of the experimental hysteresis loops in terms of local magnetization 
reversal. This model takes into account an effective uniaxial anisotropy which is the result 
of the cluster anisotropy model and a random magnetic anisotropy (RMA) as proposed by 
Harris et al. to investigate amorphous TM-RE compounds. Each Dy moment which is 
not inside a cluster, i.e. in the matrix, is assigned an uniaxial easy axis rij, these axes being 
randomly distributed. We have assumed that the anisotropy constant Dd y is the same for 
all Dy atoms. The anisotropy energy term can then be written as: 

E a = -D Dy ( £ (S,.n. t ) 2 + £ (S.z,) 2 ). (14) 

iSDy in matrix *SDy in clusters 

A. Abrupt profile 

The hysteresis loops of the multilayers with an abrupt profile (the cluster concentration 
is equal to 5%) for an in plane and normal applied field at T = IK are represented in Fig. [HJ 
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FIG. 8: Hysteresis loops of an Fe/Dy multilayer with an abrupt profile, c c i us ter =5% at T = IK for 
different Dy anisotropy constant values ((a) perpendicular; (b) parallel). The horizontal dashed 
lines correspond to M ferrimagnet ic/M sat . 
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FIG. 9: (a) Reduced magnetization of each atomic plane of a multilayer with an abrupt profile, 
Dvy/kft = 30K at T = IK for different magnetic field, (b) Corresponding schematic configurations. 
(The applied field is along the perpendicular orientation). 

For D-Dy/ k-Q = 10K and 20K, the hysteresis loops are roughly those of a system with PMA: 
the loop for a perpendicular applied field is almost square and there is no remanence when the 
field is applied in the plane (see Table [II]) • This behavior is due to the cluster anisotropy. 
The continous decrease of the magnetization before the reversal for a perpendicular field 
(Fig. [51(a)) is the result of an increase of sperimagnetism because of RMA (mainly in the 
Dy layers). The magnetization reversal mechanism is nearly uniform rotation in this case. 
For Dny/k-Q = 30K, no significant change is observed when the field is normal to the film 
except, of course, an increase of the coercive field. On the other hand, the loop is clearly 
open when the field is in the plane indicating a significant effect of the RMA. 

The reduced magnetization profile and the corresponding magnetic configurations for 
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different applied field (normal to the film) at T = IK are shown in Fig. [91(a) and Fig. 0(b), 
respectively. At B = 15T, the ferrimagnetic order is broken in the core of the Fe layer 
which leads to a magnetization larger than Mf errimagnetic as it can be seen in Fig. El (a); this 
effect has not been observed experimentally since the applied magnetic field is not strong 
enough. Note that the Fe magnetic moments in each xy plane are parallel whereas a small 
deviation from colinearity can be observed in the Dy planes (M Dy /M°[ ~ 0.95) (Fig. [3(a)). 
As the field is reduced, down to 5T, two features are responsible for the decrease of the 
magnetization in Fig. [81(a): the Fe moments tend to be antiparallel to the Dy sublattice 
magnetization and the angular distribution of the Dy moments broadens (mainly in the 
core of the layer) (Fig. [9]). The lowering of the remanent magnetization in the case of a 
perpendicular applied field compared to the cases Dvy/k B = 10K and 20K is clearly due 
to a more pronounced sperimagnetism. The dispersion of the moments is maximum just 
before the reversal (B ~ — 2T in this case). After the reversal, the angular distribution 
of the moments sharpens as the field intensity increases and the average Fe moments in 
the core of the layer deviate from the z axis. For DD y /k B = 40K, the hysteresis loops 
for a field in the plane or normal to it are quite similar showing that the magnetization 
reversal is, this time, essentially governed by the RMA and exchange interactions. This 
is consistent with the coercive field values for the two orientations reported in table. [Ill 
Indeed the coercive field for an in-plane or perpendicular applied field are very close for 
D-Dy/k-Q > 30K within the uncertainties which are significant in this case owing to strong 
competition between RMA and exchange interactions. So, our investigation evidences a 
crossover from an uniaxial anisotropy (PMA) type behavior at low Dd y values to a RMA 
type behavior at high Dj) y values. Of course, in this latter case, the magnetization reversal 
mechanism is clearly non uniform. In order to explain quantitatively this crossover, we have 
estimated the exchange and anisotropy contributions to the local field on a Dy site in the 
matrix. These two contributions are respectively given by: 

^oH^y = -^Dy-Dy-y^ X~ ^ m Dy> 



^ Dy = 7 iyDy V2 ' K Dy .n>, (16) 



2 D D y(g Dy - V 2 

(^Dy/^s) 

Assuming that the Dy moments are almost parallel to the z axis, we have HqH c 
(5ZJ Dy _ Dy )/(16/i B ) ~ 36. 5T (Z = 12) and fi H a = {5D D y) / (8fi B ) (using < m^y.n, > 
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m-Qy/2 where <> means average over the anisotropy axes). The average value of the 
anisotropy contribution varies from 9.3T to 37.2T. Actually, these two contributions are 
of the same order of magnitude when D Dy //c B ~ 30K which is consistent with the observed 
change in the magnetization reversal mechanism. 

We would like to emphasize that RMA increases significantly the coercive field values 
for large enough anisotropy constant Dd y in comparison with the cluster anisotropy model 
(Sec. UVj) . Indeed, for small values of Dv y , the magnetization reversal is roughly uniform 
and the RMA effect is very small. On the other hand, increasing Dn y leads to a non- 
uniform reversal process. Then, each magnetic moment is individually sensitive to the local 
anisotropy field which is due to RMA on the Dy sites of the matrix. 

Finally, let us mention that although the fraction of Dy atoms in the clusters relatively 
to the total number of Dy atoms is constant (3.4%), it is quite surprising that the multilayer 
behaves either as an uniaxial anisotropy system (D^ y small) or a RMA system (-Do y large). 
In the limit of very large anisotropy constant values, the small crystallized clusters would 
have no effect. 

B. APT profile 

The case of the APT profile is very interesting since it allows to investigate the influence 
of the cluster concentration. In Fig. [101 we have shown the hysteresis loops for different 
cluster concentrations at T = IK and for D^/k-Q = 30K. The main effect of the cluster 
concentration is to enhance the coercive field. The hysteresis loops in Fig. [TUJ (a) indicate an 
easy axis along the normal to the film plane. The small decrease of the magnetization before 
the reversal (for c c i uster = 5% and 10%) is attributed to the increase of the dispersion of the 
moments owing to the RMA. This effect tends to disappear as the cluster concentration is 
increased since this latter favors uniform rotation of the spins. 

When the field is applied in the film plane (Fig. [TUJ (b)), we observe that the uncertainties 
are quite large (AB C ~ 0.5T) since the coercive field does not vary monotonously with c c i ust0 r 
as it should be. Consequently, no significant effect of the cluster concentration is noted. 
However, our results indicate that the coercive field values are non zero due to the RMA but 
are much smaller than in the perpendicular orientation. As for the abrupt profile, we have 
observed that the magnetic moment dispersion is maximum just before the magnetization 



15 



-clusters 5% 
-clusters 10% 
-clusters 20% 



{Ltd** 



(a) 




B (T) B (T) 

FIG. 10: Hysteresis loops of an Fe/Dy multilayer with an APT profile at T = IK and D-Qy/k-Q 
30K for different values of the cluster concentration ((a) perpendicular; (b) parallel). 

reversal (not shown here). 
C. Discussion 



Our results evidence that the APT profile favors perpendicular anisotropy in comparison 
to the abrupt profile (see Figs. IBl and fTUi) . Concerning the abrupt profile, we observe a PMA 
behavior only when the Dy anisotropy constant value is lower than 30K. Experimentally, 
the hysteresis loops measured on the multilayers elaborated at low temperature, with a 
concentration profile close to an abrupt profile do not display any PMA unlike the sample 
built up at T = 570K. It means that our simulation model can be qualitatively compared 
to the experimental case when the anisotropy constant value D-^y/k-Q is larger than 30K. 
For Z^Dy/^B = 30K and c c i us ter = 5%, the multilayer with an APT profile exhibits roughly 
uniaxial magnetic system hysteresis loops (Bq ~ 2.00±0.15T) contrary to the abrupt profile 
multilayer. One should note that the coercive field is slightly larger for the abrupt profile 
(Bq ~ 2.45 ± 0.15T) unlike the case without RMA (see Sec. II VI) . This can be explained as 
following: as previously mentioned 24| the sperimagnetism due to RMA is more pronounced 
for the abrupt profile, so the increase of Bq when adding RMA is larger than for the APT 
profile. 
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VI. INFLUENCE OF THE TEMPERATURE 

The temperature effect has been investigated in the case of the APT profile and compared 
to the experimental results of an (Fe 3nm/Dy 2nm) multilayer elaborated at 570K. The 
simulated and experimental hysteresis loops at T = 5, 50 and 100K are shown in Fig. [TTJ 
The numerical results have been obtained with D-Dy/ks = 40K and c c i ustC r =10% for which 
the best agreement with experimental data at low temperature has been observed. At 
T = 5K, the multilayer exhibits PMA unlike the multilayer with an abrupt profile (previous 
section). Note that the loop with an in-plane field at T = 5K is open contrary to the case 
D-Dy/k-Q = 30K (previous section). 

The temperature dependence of the simulated loops is in good qualitative agreement 
with the experimental ones. For each temperature, the hysteresis loops evidence an easy 
axis normal to the film. The coercive field decreases as the temperature increases because of 
the lowering of the effective anisotropy constant in relation with thermal fluctuations of the 
atomic moments. The hysteresis loops with an applied field in the film plane are open at low 
temperature owing to RMA. Above 50K, no hysteretic phenomenon is observed indicating 
that RMA has no effect anymore in this temperature range. The hysteresis loops of our 
model above 50K are those of a uniaxial anisotropy type system with an easy orientation 
along the z axis. It can be seen that unlike the simulated loops, the experimental ones 
do not saturate mainly for an applied field in the film plane. This indicates that the fan 
structure of the Dy moments is more rigid in real systems than in our model. 

VII. CONCLUSION 

In this work, we have investigated a model of amorphous Fe/Dy multilayers with RMA 
and small crystallized clusters defining on average a preferential axis normal to the film plane. 
Our results using a model built up from experimental concentration profiles are consistent 
with measured hysteresis loops, i.e. APT profile favors PMA in comparison to the abrupt 
profile. It has been shown that the crystallized clusters embedded in the matrix with a 
very small concentration could explain the observed PMA. Moreover, the cluster anisotropy 
effect is more pronounced with the APT profile. With the abrupt profile, RMA plays a more 
significant role and may even hide the cluster anisotropy effect when the anisotropy constant 
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FIG. 11: Hysteresis loops of Fe/Dy multilayers at T = 5,50, 100K. (a, c, e) simulated loops with 
an APT profile, ^Dy/te = 40K and c c i uster =10%; (b, d, f) experimental loops for multilayer 
elaborated at 570K. 

is enlarged. The temperature influence on a multilayer with an APT profile is to remove 
RMA effect. In a near future, we plan to investigate magnetostrictive and magnetoresistive 
multilayers. 
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TABLE I: Coercive field values at T = IK of the multilayer with the abrupt and the APT profiles 
with Cduster = 5% for D-Dy/k-Q = 10K and 50K. 



D Dy /k B (K) 


B c (T) Abrupt 


B c (T) APT 


10 


0.41 ±0.15 


0.57 ±0.15 


50 


1.38 ±0.15 


2.13 ±0.15 
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TABLE II: Coercive field values at T = IK of the multilayer with the abrupt profile with c c i ustC r = 
5% for various values of DDy/^B hi perpendicular and parallel field orientations. 



^(K) 


10 


20 


30 


40 


Bh (T) 


0.33 ±0.15 


1.10 ±0.15 


2.45 ± 0.30 


9.80 ±0.50 


b'A (t) 


0.10 ±0.15 


0.20 ±0.15 


2.20 ±0.30 


8.40 ± 0.50 
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